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We continue the investigation of massive integrable models by means of the bootstrap fusion 
procedure, started in our previous work on 0(3) nonlinear sigma model. Using the analogy with 
SU(2) Thirring model and the 0(3) nonlinear sigma model we prove a similar relation between 
' sine-Gordon theory and a one-parameter deformation of the 0(3) sigma model, the sausage model. 

^\ , This allows us to write down a free field representation for the Zamolodchikov-Faddeev algebra of 

0^ ' the sausage model and to construct an integral representation for the generating functions of form- 

factors in this theory. We also clear up the origin of the singularities in the bootstrap construction 
(— I . and the reason for the problem with the kinematical poles. 
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I. INTRODUCTION 

Exactly solvable two-dimensional quantum field theories have attracted a great deal of interest and research work 
^ , recently. There are several reasons for that: these theories can serve as laboratory for higher dimensional theories; 
• there are physical systems with effectively two-dimensional phenomena (eg. quantum Hall effect); and there is the 
huge and rapidly developing area of string theory, topological field theory and related topics, in which two-dimensional 
physics plays a central role. 

A large variety of different approaches emerged in trying to solve these models of quantum field theory. The first 
real breakthrough was achieved in in which exact S-matrices were conjectured for several interesting models using 
bootstrap ideas. The fact that made this possible is that scattering theory of quantum integrable models is very 
simple: it is factorized scattering theory (FST). This means that the many-particle S-matrix can be written in terms 
of products of two-particle S-matrices. The consistency condition for this turns out to be the famous Yang-Baxter 
equation, which, together with unitarity, crossing symmetry and the spectrum of the particles of the model, severely 
restricts the class of possible S-matrices that one can write down. However, there is a so-called CDD ambiguity, which 
means that the S-matrix is not determined uniquely by the aforementioned conditions. In most cases one solves this 
difficulty by choosing a 'minimal' solution to the factorization equations. This choice can be justified in many cases 
by comparison of results coming from the conjectured S-matrix and from the Lagrangian field theory. 

The S-matrix, however, describes only the on-shell physics. For off-shell physics, one turns to other quantities 
5J] ' which contain more information: correlation functions. For their calculation the knowledge of the matrix elements of 
I local operator is needed. These matrix elements are called form-factors; correlation functions can be represented as 
an infinite series constructed using these objects. For the majority of integrable field theories this is the only tool to 
compute correlation functions up to now; in some simple cases one can use other techniques. 

Form-factors can be found solving a very complicated system of axioms, which have been known for a long time 
The two-particle minimal form-factors were calculated for many models, however the calculation of the many-particle 
ones posed essential difficulties. Then, in a remarkable series of papers, Smirnov obtained all form-factors of 

some important operators in the sine-Gordon, SU(2) Thirring and 0(3) sigma models. Further investigation revealed 
that the key structure behind these solutions is that of quantum affine algebras, which is a nonlocal symmetry of 
these integrable field theories 1^. The Zamolodchikov-Faddeev operators (which play a central role in our approach 
as well) turned out to be vertex operators for the representations of the nonlocal symmetry algebra ^ . 

The analogy between the Knizhnik-Zamolodchikov (KZ) equations for affine Kac-Moody algebras and the form- 
factor axioms, one of which proved to be a quantum analogue of the KZ equation HJl^, as well as the vertex operator 
representation of quantum affine algebras |14|| strongly suggest that a bosonization technique based on free fields, 
similar to the one used in conformal field theory, should be useful in solving the form-factor equations. This is the 
central idea in the papers of Lukyanov [p^|jl6| , where he was able to construct a free field representation for the sine- 
Gordon and SU(2) Thirring model and to obtain an integral representation for the form factors. In our earlier work [Q, 
using mainly his results and the connection between SU(2) Thirring and 0(3) sigma models put forward by Smirnov, 
we defined a free field representation for the 0(3) model and applied it to build up an integral representation for the 
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form- factors. Here we extend the method to the sine-Gordon and sausage models and in this way write down a similar 
representation for the form-factors of the sausage model. The problems encountered in |l) when applying the method 
unfortunately remain the same and we choose the same old solution for correcting the analytic structure introducing 
some multiplier functions. Our main result is to obtain a procedure following which one can calculate form- factors of 
the sausage model. We also provide a detailed analysis of the analytical structure and an interpretation of the double 
poles in terms of an analogy with thermal Green functions. The relation of the problem of the kincmatical poles and 
the singularity of the fusion limit of the sine-Gordon form-factor is cleared up. 

The paper is organized as follows. Section II reviews factorized scattering theory of sine-Gordon and sausage models 
and there we discuss the quantum symmetries. In Section III we briefly describe the form-factor axioms and the free 
field representation for their solution. Section IV deals with the question of the bootstrap fusion procedure used for 
transition from the sine-Gordon model to the sausage model and using these results defines the free field representation 
for the sausage model, while in Section V we give the general formulas for evaluating the form-factors and provide a 
detailed discussion of their analytic structure. Section VI is reserved for the conclusions. 

We added three appendices to the paper. Appendix A describes the method used to obtain the bootstrap fusion 
and the Zamolodchikov-Faddeev algebra relation on the level of the free field representation, while Appendix B shows 
how to evaulate the traces occurring in the calculation of generating functions with coherent oscillator methods. In 
Appendix C we describe a sample calculation of the operator product singularities. 

II. REVIEW OF THE SINE-GORDON AND SAUSAGE FIELD THEORIES 

Sine-Gordon theory is defined by the Lagrangian density 

C^^id.ct^f + '^cosihcj^), (2.1) 
where 6 is a coupling constant and toq is a mass parameter. Instead of b we will use the parameter 

which can be taken as the only dimcnsionless parameter of the model. When ^ > 1, the quantum spectrum of the 
theory contains only a doublet of solitons which forms the fundamental representation of the quantum group SU{2)q 
with 

q = exp[iTT{l + i)] . (2.3) 

When ^ < 1, the spectrum contains breathers, which are bound states of the fundamental solitons. We do not treat 
this case here. We also restrict ourselves to study the case when q is not a root of unity to avoid dealing with 
degenerate representations of the quantum group. The S-matrix of the solitons was found in ||^ and looks as follows: 



S+tiP) ^ SZZ{(3) = SoW), 



sinh(f^ 
smh( ,^ ) 



sinh(2|) 
sinh(H^) 



Sl:^iP) = S^li(3) ^ SoiP) , IJi (2.4) 



where — is the internal index, (3 is the rapidity parameter, which is related to the energy-momentum by pq 
mcosh/3, pi = msinh/3, m is the soliton mass and 



•S'o(/3) = -exp 







°° i-f 

sin(K/3) smh{^^nK)dK, 



K cosh(i7rK) sinh(|7rK) 



(2.5) 



The Hilbert space of the theory can be described by introducing the formal Zamolodchikov-Faddeev operators Zi{P) 
satisfying the following algebra relations 
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Z,{f3i)Z02) = S^^iPi2)ZiiP2)Zk{f3i), /3i2 =Pi- 132. (2.6) 

As we shall see below, t he c ommutation relations of the Z-operators reflect the scattering of the particles. From 
now on we will refer to ( |2.6| ) as the scattering relations for the sine-Gordon model. The Hilbert space furnishes a 
representation of the Zamolodchikov-Faddeev (ZF) algebra. The ZF operators create the asymptotic in and out states: 

\ZM)...Z,^{(3n))out = Z,,(/3i)...Z,„(/?„)|vac), 

|Z,,(/3i)...Z,„ (/?„))„ = Z,„(/3„)...Z,,(/3i)|vac) , (2.7) 

where |vac) is the vacuum state and the rapidities are ordered as /?i < /32 < •■■ < Now if one expresses the 
out states in terms of the in states, one gets the n-particle S-matrix as a product of the 2-particle S-matrices. The 
associativity of the ZF algebra then implies the famous Yang-Baxter equation, which is the factorization condition 
for the S-matrix ||]. 

The Hilbert space structure is specified by defining the conjugate ZF operators Z'^\ which annihilate the vacuum 
and satisfy the following additional relations 

Z^^iP,)Z^\f32) = Z'H(32)Z'^HP^)SI]{(3,2) , 

Z,(/3i)ZJ"t(/32) = Z''Hp2)Sg{f3i2)Z,{f3i) + 27rSfS{(3i2)- (2.8) 



However, in the sequel we will deal with another representation of the ZF algebra defined by (2^) which does not 
include the conjugate ZF operators and therefore it is different from the physical Hilbert space of the model. This 
representation will be the space of the free field system used to solve the form- factor equations. 

As mentioned above, the sine-Gordon model is invariant under SU{2)q quantum group symmetry. To make this 
symmetry manifest, one has to redefine the ZF operators 

Za{P)^exp(a^)Za{f3) , (2.9) 



2^ 

which implies the following redefinition of the S-matrix 

S:t{P) = exp (^{a c)QsiiiP). (2.10) 
It can be easily proved that if we identify the vectors Za(/3)|vac) as the fundamental representation of the quantum 



group SU{2)q, then the redefined S-matrix and the scattering relations (2.6) are invariant under the action of the 
quantum group. However, care must be taken since S is not a unitary matrix and the representation of the out states 
differs from that of the in states by a change q — > 1/q. Since \q\ = 1, this is equivalent to a complex conjugation. Using 
this invariance and the fact that the decomposition of the corresponding representations is given hy^®^—0®l (we 
label the representations using the conventional isospin notation), we can extract the amplitudes in the two possible 
channels which turn out to be 



Sisospin-l{P) = 'S'o(/3), Sisospin-o{P) = -p _i_ /3\ ^oiP) , (2-11) 

sinh ( IZL+iZ 



and S can be decomposed as 



sinh(^) 

T 

S — Sisospin — O-Po Sisospin — l-Pl ; (2T2) 



where Pq and Pi are projectors on the spin-0 and spin-1 subspace of the product of two spin- 1/2 representations. 

The sausage model was introduced in JlTf as a one-parameter family of deformations of the 0(3) nonlinear sigma 
model. The spectrum of the model consists of a massive (quantum group) isotriplet of scalars, capital letters /, J, K, L 
will be used to label the internal quantum number which takes the values -l-,0, — . The S-matrix is invariant under a 
U{1) symmetry and the charge conjugation which takes I = —I. The S-matrix has the following symmetries: 

Sfj^iP) = Sfj\p) = (/3) = SlidP) (C,P,T) , 

5f/(/3) = S^jf{i^ - 13) (crossing) . (2.13) 
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The independent S-matrix-elements are 



sinh A(/3 — in) 



sttW 

sin 27rA 



sinh A (/3 + in) 



sin ttA sin 2nX 



sinh A(/3 + itt) sinh A(/3 — 2in) 
sinh A/3 
sinh A (/3 - 2in) ' 



5o"°(/3)=5+°(/3)+5+l(/?) , (2.14) 

where A is a parameter of the model. This is a consistent S-matrix if < A < 1/2; otherwise, there are poles inside 
the physical sheet < Im/3 < n. In the one-loop renormalized action can be found. For our considerations we 



will need only the S-matrix ( 2.14 ) 



The S-matrix of the sausage-model is SU{2)q invariant with q = exp[i7r(A + 1)]. When we set A — > 0, we find 
the S-matrix of the 0(3) nonlinear sigma model. One can introduce the ZF operators 2j{f3) for the sausage model 
similarly to the case of sine-Gordon theory. The basis needed to exhibit the quantum group symmetry is given by 

i/(/3) =exp(A//3)Z/(/3). (2.15) 

The states created by the operators Z form a triplet (isospin-1) representation of SU{2)q. Similarly to the case of 
sine-Gordon model (see ( p. 12 ) ) , we can decompose the S-matrix into three channels corresponding to 1®\ = 0©1©2. 



The amplitudes for the three channels are 

5o(/3) = 



sinhA(/3 + 2in) 

sinhA(/3- 2i7r) ' 

sinh \{(3 + 2in) sinh A(/3 - in) 

sinh A(/3 — 2in) sinh X{P + in) 



sinh X{P- in) 

^'^^^ ^ sinhA(/3 + z.) • ^'-^^^ 



III. FORM-FACTOR AXIOMS AND THE FREE FIELD REPRESENTATION FOR THEIR SOLUTION 



In Smirnov's approach to integrable field theories, the local operators are defined in terms of their form-factors. 
Form-factors are the matrix elements of the local operators between the vacuum and many-particle states: 

F^°...,^JPu■■■,Pn) = (z^acj 0(0) ,...,/?„),„.. . (3.1) 

By virtue of crossing symmetry, these matrix elements determine the general matrix element between two many- 
particle states and, assuming asynptotic completeness of Hilbert space, this means the complete characterization of 
the local operator 0{x). The form-factors of a local operator must satisfy the following set of axioms: 

1. The functions F are analytic in the rapidity differences /3y = (3i — (3j inside the strip < Im/3 < 2n except 
for simple poles. The physical matrix elements are the values of the -F's at real rapidity arguments ordered as 
/3i < /32 < . . . < 

2. The form- factor of a local operator of Lorentz spin s satisfies 

F^,,...,^SPl A, ...,/?„ + A) = cxp(sA)F,,,...,,„(/3i, . . . , /3„). (3.2) 

This means that the form-factor of a spin s — Q operator depends only on the rapidity differences. 

3. Form-factors satisfy Watson's symmetry property: 

Fii,...,ij + iPl : ■ ■ ■ I /3j , Pj+l j • ■ • j Pn) — 

S-^t^ F,,_k,^,,u„...,^APl, ■ ■ • ■ • (3.3) 
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4. Form-factors satisfy the cyclic property 

(/3i, ...,/?„ + 2TTi) = (/?„, . . . , /?„_i). (3.4) 

In the case of an operator which has a mutual locality index (jj(0,5') with respect to the elementary field 5*, this 
axiom is modified to 

F^„...,^„ {pi,...,(3n + 2«) = exp[2^iw(0, (/?„, . . . , (3n-l). (3.5) 

5. Form-factors have kinematical singularities which are simple poles at points where two of the rapidity arguments 
are displaced by in. The residue is given by the following formula: 

res/3_^„+i^Fji^...^i„(^i, . . . ,/5„) = -. (^i, /3j, ...,/?„) x 

J 17 - ■ ■ J j " "> n — l 

We would like to remark here that the cyclicity axiom closely resembles the KMS condition for finite-temperature 
Green functions of quantum field theory. We will see that the expressions for the solutions of the form-factor equations 
really look like expectation values in a heat bath and later we will use this analogy. 

The locality theorem proved in ||^ ensures that the operators defined by solutions of the above equations are local 
(up to the mutual locality indices, if there are such). The importance of the form- factors is given by the fact that 
the correlation function of two local operators can be expanded into an infinite scries in terms of their form-factors, 
obtained by inserting a complete set of many-particle states into the correlator. 

In the paper of Lukyanov , a free field representation for the Zamolodchikov-Faddeev algebra was proposed to 
solve the above axioms. That work considered the SU(2) Thirring and the sine-Gordon models. Here we recall some 
essential facts about this free field realization. We do not give a full treatment and restrict ourselves to the formulas, 
used in describing sine- Gordon form-factors, which are necessary in the subsequent sections. The interested reader is 
referred to [0 for more details. 

Let us introduce a field <j){(3) with the following properties: 

[m),</'(/32)] =ln^o(/32-/3i), 

(O|0(/3i)0(/32)|O) = -ln5(/52-/3i), (3.7) 



where iSo is the function defined in eqn. ( p.5[ ) and we give the formula for the function g below ( p.ll| ) . The consistency 
of the two relations above requires 

The field (f> is represented (after a proper ultraviolet cut-off procedure) on a Fock space. We will build up operators 
which represent the ZF algebra on that space. 
The field 

^(/3)-(/)(/3 + z|) + 0(/?-z|), (3.9) 
will be important to define the screening charge below. Its two-point functions can be written down as 

(O|^(/3i)0(/32)|O) =lnu;(/32-/3i), 

(0|^(/3i)<^(/32)|0) = - \ng{P2 - Pi), (3.10) 
where we defined the following functions 



9^P} ^r(i)^ r(f) MP) 



, r(-^ + f) 

r(i + ^ + f) ' 
^r(i+|+^) 
< r(-i + f) 



w(/3) = K 
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where k is a normalization constant and Ri is the function 



Now we define Lukyanov's vertex operators: 

Vif3) = exp(z0(/3)) - ig{0))i : exp(i(/.(/3)) :, 

= exp(-i0(/3)) = (g(0))^ : exp(-i^(/3)) :, (3.f3) 

where : denotes an appropriate normal ordering [|6| . These operators need to be regularized because the functions g 
and g have simple zeros at f3 — 0. Using the analogy with conformal field theory we define the regularized values of 
the functions g,g at /? = to be 

r 2 

9reg[0) = ^™ — =' P ' 

5.e,(0)= lim^=:p-^ (3.14) 

The vertex operators defined in ( ^.13| ) obey the following fundamental relations: 

V{P^)V{p2) = P'9{P2 - Pi) : V{Pi)V{p2) : , 
V{(3i)V{(32) = pMP2 - Pi) ■■ V{f3i)V{f32) : , 

V{f3i)V{f32)=p^g{p2-pi):V{f3i)V{f32): . (3.15) 

Eqn. ( 3.15| ) plays a key role in proving the ZF scattering relations. 

Then in analogy with the Coulomb gas representation of rational conformal field theory one introduces a screening 
charge x by the definition 

{u\x\v) ^ V-\u\^ f d^V{^)\v), (3.16) 

where C is a contour specified in the following manner: assuming that all matrix elements (u|y(7)|w) are meromorphic 
functions decreasing at the infinity faster than 7^^, the contour goes from 5ft7 = — 00 to 3?7 — +00 and lies above 
all singularities whose positions depend on \u) and below all singularities depending on \v). {rj is a normalization 
parameter which is chosen to get the correct residue for the pole on the operator product below (see [l^)). 
Having defined these objects, one can write down the following operators: 

Z+(/3)-exp(-^)F(/3), 

Z_(/?) = eM^)i<lhv{P) - q--^ViP)x). (3.17) 
Then the following propositions hold |l6| : 

(i) These operators satisfy the SU(2) Thirring model ZF relations. 

(ii) The singular part of the operator product Zi{P2)Zj{f3i) considered as a function of the complex variable [^2 for 
real /3i in the upper half plane 3/32 > contains only one simple pole at (32 = Pi + ™ with residue —iCij. 

We introduce the following structures, specified later in details: 

1. The operator of Lorentz boosts K satisfying 

Z,;(/3 + A) = exp(-Ai^)Z,;(/3) exp(Ai^). (3.18) 

2. A map O L{0) from the space of local operators to the algebra of endomorphisms of ttz satisfying the 
following two conditions 

L{0)Z,{f3) = Z,{P)L{0), 
exp{AK)L{0) cxp{-AK) ^ exp{As{0))L{0). (3.19) 
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where s{0) denotes the Lorentz spin of the operator O. 

If one has these objects, then, as shown by Lukyanov |J6|, the functions 

. . . ,/3„) = rr,^[exp(27rz/0L(O)Z,J/?„) . . . Z^)], (3.20) 

are sohitions of the form- factor axioms. 

We now have to give the definition of the dual operators which give us the map L representing the local operators 
of the model. We introduce the scalar field satisfying the same relations as 4> together with its associated field 

(j)'{a). (We give a closer definition of these fields shortly when we introduce an ultraviolet regularization for the free 
field system.) With the help of these new fields we define the dual fields as follows: 



A_(a) = U'{a), 



Ao{a) = - Wx'U'ia) - q'-^U' {a)x'], 

^2cos(^) 

^+(") = 9rJ( ^ ^ il'x''u'ia) (q' + q'-')x'U'{a)x' + q'-'U'ia)x'% (3.21) 

Z COS^^qTj-J 

The dual vertex operators U' and the V' required for satisfy similar relations as the unprimed ones, but with the 
function g changed to 

-(^ + 1) r((^ + ^) 

The parameter q' equals exp 

One can also cj 
using the relation 



One can also calculate all the normal ordering relations involving the primed and unprimed operators merely by 



V'{d)V{P) = pp'uiP - S) : V'{d)V{P) :, (3.23) 

where the function u is given by 

u{f3) = k'^, (3.24) 

with k a constant. 

Now one writes down the functions 

J''j\'::2{ai, . ..ak\Pi, . . ./3„) = Tr^^[exp(2^^i^)A,,(afc) . . . A,, (ai)Z,, (A) . ..Z,,{l3i)], (3.25) 
then these functions have to be expanded as 

(ai,. . .afc|/3i, . . ./3„) = ^ F'^^-^-iai, . ■■au\{sj})Fj,_j,{{sj}\l3i, .../?„) x exp(siai) . . .exp(sfcafe), (3.26) 

where the functions F are the required solutions of the form-factor axioms. This means that the functions T play the 
role of generating functions for the form- factors of local operators. They describe n-particle form-factors of infinitely 
many operators which may be considered as descendents of the first one in the sense of j|] . 

The actual calculation of the matrix elements and the traces proceeds through the procedure described in . The 
free field construction needs a regularization; we choose an ultraviolet cut-off by taking the rapidity interval finite 

TT TT 

--</?<-. (3.27) 
e e 



Then we introduce the mode expansion of the free field (p as: 



m = V~siQ - iPP) + ■ ■ ■ f . ^ exp(ifce/3), (3.28) 
^ — ' I smh TTKe 



with 
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s = (3.29) 



The oscillator modes satisfy the commutation relation 

sinh 2^ sinhTrfce smh^(£ + 1) 

I.,.,]. — 1_ -f4i^**.- ("0) 

while for the zero modes we have the usual canonical ones: 

[P,Q]=-i. (3.31) 

The two-point function and the commutation relations for this field can be calculated and it can be seen that in the 
limit i we recover the relations satisfyed by the field (see eq. (3.7)). The dual field is given by 

<^^(a) = -V7(g - laP) + J2 ^-^frr ^M^kea), s' = ,^7^y (3.32) 

^ismh(7rfce) 2(^ + 1) 

with the oscillators a'^ defined by the relation 

a^, sinh— (^ + 1) = Ofcsinh— ^ . (3.33) 
We introduce the conventional notation used in conformal field theory 



«+=-l/a_ = ^^. (3.34) 

The dual set of operators can actually be obtained by the substitution S, ^ —I ~ ^, which is known to have an analogy 
in conformal field theory (see Lukyanov for details) . The space where the operators and their duals act is defined 
to be 

ttI - , (3.35) 



V2 



where Fp denotes the Fock space built up with the help of the creation operators (as usual, the oscillator modes with 
negativ indices) from the ground state \p) which satisfies: 

P\p)^p\p). (3.36) 

The operator of the Lorentz boost is given by 

K, = leH + -±-—P^ (3.37) 
V2 

where 

+ 7 J—- a'-fcOfe. (3.38) 

2 f — ' sinh 2^ sinh irke 

k—l ^ 

Formally we can write tt^ — lim tt^. This should be understood in the following sense: we make all calculations with 
the regularized operators and in the end take the limit e — s- 0. 



IV. SAUSAGE MODEL FROM BOOTSTRAP FUSION OF SINE-GORDON MODEL 

In this section we describe how to obtain the sausage model S-matrix introduced in section || from that of the 
sine-Gordon theory. First we introduce some graphical notations. We denote the S-matrix by a cross: 
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This should be understood as foUows: the hues are vector spaces carrying the internal indices of the particles. The 
S-matrix acts in the tensor product of two such spaces and carries an argument which is the rapidity difference of 
the two scattering particle. The rapidities are associated to the lines of the diagram. Each segment carries an index 
and repeated indices must be summed over. The order of the matrix products is represented by the arrows. As an 
example, the Yang-Baxter equation in these notations can be written as (omitting the rapidities of the lines, which 
causes no confusion) 





which means 

4f (/323)5g(/5i3)^:A:|(/3i2) = 5g(/3i2)5:^f (/5i3)5:i:|(/323). 

The bootstrap fusion means the following formal operation on the level of the ZF operators 



(4.1) 



(4.2) 



where (7)^*^ denotes the appropriate Clebsh-Gordan coefficients for the internal indices. In our case they are just the 
coefficients for the spin-1 representation in the product of two spin-1/2 representations. The nonzero elements are: 



C++ = cz- = i, c+- = 



1 



-1/2 (J- 



7I/2 



(4.3) 



Eqn. (4.2) is formal since the rapidity difference just corresponds to the position of the annihilation poles and so there 
are potential divergences. Fortunately, in the case we investigate thes e singularities drop out and we will be able to 
define the ZF operators corresponding to the RHS of the equation |4.2| in the free field representation. We will denote 
this operation graphically by 



> 



I 



Then we can write down the graphical representation for the S-matrix obtained from the bootstrap fusion: 




'5///|(/3i2) 
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This diagram s hou ld be read as follows: two incoming particle in the triplet representation is converted (using the 
fusion formula into two pairs of doublets and then scattered on each other (this means four crossings ie. four 
S-matrix factors). After that, the triplet part is projected out using again the fusion formula. However, care must be 
taken, since in the case of the incoming particles we make an embedding of the triplet representation into the product 
of the two doublet representations, while in the case of the outgoing particles we go the other way via projection. 
These two operations are adjoint to each other and differ not only by transposition of matrices, but (and here comes 
an essential difference from the 0(3) sigma model case treated in |l|) also differ by a substitution q ^ 1/q, which 
is here just a complex conjugation since |g| = 1. We can write the fusion equation in the following form (using the 
redefinitions 2^, 2.10| and |2.15 ) 



(4.4) 



The coefficients C and C are given in (4.2). We remark that it is necessary to use the redefined S-matrices because 
to apply the quantum group Clebsh-Gordan coefficients we have to go to the one-particle basis which is at the same 
time the properly normalized basis of the quantum group representation. 

The proof of the fusion relation between the sine-Gordon and sausage model S-matrices proceeds as follows. The 
sine-Gordon S-matrix acts in the tensor product of two spaces; we can expand it in terms of tensor products of Pauli 
matrices acting in one space as follows: 



1 



2 sinh( 



(3 — ITT > 



/3- 



sinh 



) sinh 



sinh ( ^ 
■ s 



> a\ ^ sinh (^~~7~~) + {~? 



(4.5) 



The redefinition means changing the Pauli matrices by rapidity-dependent factors. The product of the four S-matrices 
can be written in the tensor product of four two-dimensional spaces, each S-matrix acting only in two of them: eg. 
the first S-matrix in the spaces 1 and 4, the second one in the spaces 1 and 3 and so on. Then one can embed 
the Clebsh-Gordan coefficients into this 16 dimensional space and carry out the matrix multiplication. The special 
structure reflected in the graphical presentation of the formula greatly simplifies the required computation. The 
relation between the parameters of the two models turns out to be A = 1/^. 

Now we turn to the fusion of the ZF operators. We want to give a meaning to the formula |4.2| . We interpret it as 
the limit 



Zj(/3) = lim Cr^Z,, (/?! + -+ 5)Z,, {(32 
o — *o z 



ITT 

y 



<5). 



(4.6) 



Remarkably, all the singularities cancel and the result is a very simple formula for in terms of the free field represen- 
tation: 



exp 



2cos(f ) 



cos(f)^^nC 



(I) kx'UiP) - (g + q-')xUif3)x + q-'U{f3)x' 



(4.7) 



The proof of this result can be found in Appendix |^, where we also outline the direct calculation of the scattering 
relations. They prove to be identical to the sausage model ZF algebra: 



Z,(/3i)Zj(/32) = 5f/(/3i2)Zi(/32)ZK(/3i). 



(4.8) 



However, just as in the case of the 0(3) sigma model [Q, the operator product does not have the right kinematical 
poles. Instead of having them at rapidity difference itt, it turns out that the singularities are situated at 2i7r (see 
Appendix ^). Nevertheless, if we write down the trace formula 3.25, it will solve the cyclicity (by the definition of 
the trace) and Watson's axiom. The analytic structure will be corrected by a procedure similar to that in [||. 
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V. EVALUATION OF THE EXPRESSIONS FOR THE FORM-FACTORS 



We have to calculate the traces required for obtainig the integrand of the representation for the form- factors. Using 
the coherent state representation for the traces, we find the result: 



Tr^, [exp(27riif )C/'(afc) . . . U'{ai)V'{5p) . . . y'(5i)[/(/3„) . . . t/(/3i)y(7,) . . . ^(71)] 



i — np 2 p 2 ^2 



I T] 7] ^ OnrOpk exp | - 



E'^j~E"j 



n—nn 2 n 

n - n - n ^•^'(7^ - n ^'(^^ - n ^'(-^^ - '^^■) n ^"'(^^ - "^■) 
n - «,) n - -5,) n ^(7. - s,) n ^-'(7. - 



(5.1) 



^,3 



Here we defined the new functions as: 



G(/?) = -^^sinh(^±^)sinh/3, 



G'(a) = -C^ 



i/(a) 



sinh ( 



(e + l)sinh(^)^ 



and the constants are given by: 

Ci = exp 

C2 = exp 

C'2 = exp 



cosh a ' 



t 



(5.2) 





r 00 



^^sinh ^exp(-t)g-^^^^(^_^y 
t sinh 2i cosh t sinh 

2 t 

^^sinh ^exp(-t)g.^^^(^_^) 
t sinh 2t cosh t sinh 



sinh^ ^ exp(-^) g-^t^^^ 
t sinh 2t cosh t sinht(^ + l) 



(5.3) 



The trace over the zero modes contains an infinite constant which we have set equal to unity. 

We also have to specify the integration over the variables of the screening charges. We face the following problem: 
when we take the integration contour as given in the sine-Gordon model (see in [|l6| ) and try to take the limit which 
gives the sausage model, new double poles arise in the integration over the % variables whose positions depend on 
the Pj variables. This phenomenon was already noticed by Smirnov Q for the case of the 0(3) sigma model. The 
contour gets pinched by the two poles approaching each other and the integrals diverge. This means that the limit 
to the sausage model is singular. To deal with the integral we introduce the prescription that one should take the 
coefficient of the most divergent term in the integral. 

We are now to give a justification of the prescription and explain the appearance of double poles in the "finite 
temperature" expressions, while they are not present in the vacuum expectation values, ie. the "zero-temperature" 
matrix elements. 

The basic picture is as follows. Taking the trace means imposing the periodicity equation (or KMS condition) 
included in the form-factor axioms. What happens is that while in calculating matrix elements we use Green functions 
of the scalar field on the plane, here we impose the periodcity condition. The original function g{f3) has poles at the 
rapidity difference and iir, which will have all their images under the KMS and the Watson's symmetry condition. 
They can be viewed as mirror charges introduced to satisfy the boundary conditions for G(/3), which is the analogue of 
finite-temperature Green function, if we take g to be the analogue of zero temperature one. Let us for simplicity take 
the case of the 0(3) model, which means ^ = oo. We deform the trace by substituting exp(27riif) by exp{{2TT + ()iK) , 
where ^ is a small positive number. This means a deformation of form-factors which is the same as considered in [ p^ . 
The deformed form-factor wiil obey Watson's symmetry property and a deformed KMS condition: 
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...,/?„ + (27r + C)z) = J^.„,n,.. ■ ■ • , /3„-i). (5.4) 
Looking at the explicit form of the function G, when C > 0, which is (up to a normahzation factor) 



i/3 + iTr) smh(^^ 



n/3 



Y (i _ IL^llll] Y (i - IL±li 



T 

TT - 

271 + " V" 2FTC 

(for the calculation see Appendix |^), we see that the double pole in G^^ (which is relevant for the integrand of the 
contour integral over 7^) arises when ^ — > 0. However, the operator product singularities are the same for finite ^, ie 
we have a pole at P12 = 27ri and all its periodic images, and the generating function is regular at P12 — T:i. So we 
take the leading singular term of it when ( ^ 0. It will automatically obey the (undeformed) KMS condition and 
Watson's property. But its behaviour at the positions of rapidity difference in changes. In Appendix ^ it is shown 
that the nonsingular behaviour of the operator product at these points is to be attributed to a zero in the Green 
function cancelling the pole coming from the pinch of the contour integral. But the same effect, which doubles the 
pole in the integrand, also doubles this zero in the factors G{Pi — f3j) and so the net result is a function having simple 
poles at Pij = 0, 27ri and simple zeros at (3ij = iri. This is displayed well in the structure of the generating function 
for the 0(3) current, calculated in ||l[, which we remind here: 



2i7r(/92i + «7r) 
(/321 + 2z7r)/32i 



1 



cosh(/3i — a) cosh(/32 — a) 



tanh^ ^ . (5.6) 



To deal with this problem, we saw in [Q, that we had to introduce a 2ni periodic, symmetric function, which killed 
the poles at and 27rj and turned the zero at iri into a pole. This was the function (here written for n-particle case) 

n tanh^^. (5.7) 

l<i<j<n 

The end result for the generating function of the two-particle form-factors of the current (and its descendents) is 
therefore 



8 ^ /3i2(2i7r-/?i2) ^ T Vcosh(a-/3i) cosh(a - P2) j ' ' 

taking into account the proper normalization factors to get the correct value of the charges on particle states. In the 
light of the argument given above, we expect this to be a general prescription to get the right structure of kinematical 
poles. 

Although in the above discussion we took ^ — 00, it carries over to finite values of ^. In this case additional poles 
are introduced (both in g and G) but they are outside the strip since we choose ^ > 2 in order to be in the regime 
where the sausage model has no bound state poles. Taking suitably small C, it is possible to draw the contours in a 
way which avoids these additional singularities. 

For the contours of integration over the variables 5j , no similar problem arises and we have to draw the contours just 
as for the vacuum matrix elements but they must lie within the strip —in < Sj < in. This completes our prescription 
of how to write down generating functions for the form-factors of the sausage model. 

We would like to note that in the case of finite ^, the correct analogue of the above generating function is not just 
a simple continuation to finite ^ |^]. The expression 

Tr,,[A_(a)Zo(/?2)2+(/3i)] (5.9) 

turns out to be zero due to the Kronecker delta arising from the zero mode trace. The correct analogue would be to 
take eg. 

Tr^AMa)2+{P2)Z-Wi)]- (5.10) 

If one writes down this expression, one gets a complicated integral formula which however can be seen to yield the 
correct 0(3 ) va lue in the limit ^ ^ (as well as all other formulas for form-factors written down using the general 
expression (^.ij)). 

The case ^ = cxd is special since in that case the appropriate factor is S2n-2r+2p-2k, which allows much more 
nonzero solutions lllll. The reason for that is that the zero-mode trace enforces momentum conservation for the free 
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field. While in the 0(3) case a+ = — a_ and the momentum cancellation is possible between the vertex operator and 
their duals, for generic ^ this is not the case: they have to cancel separately. Since we have actually Two quantum 
groups, one with q = exp iiTa\ for the vertex operators Z and one with q' = exp i'Ka^_ for their duals, the condition 
which is imposed is actually that we should have operators under the trace which are singlets separately under each 
of them. This is very much in the spirit of Smirnov's work j^] and resembles his notion of the "main set of solutions" 
(which are the one which have their counterparts when deforming away from q = —1). If ^ is rational, however, some 
cancellations happen and the space of the solutions increases, suggesting an analogy with the reduction of sine-Gordon 
theory at rational values of the coupling, corresponding to some well-known changes in the representation theory of 
the quantum group when the value of q equals a root of unity pO| . 



In this paper we were able to show that there exists a free field representation for the sausage model, very similar to 
the one found in for the 0(3) case. We could also identify the origin of the problem with the kinematical singularities 
and the singularity in the fusion limit: the first is related to the fact that the operator product singularities of the 
fused ZF operators are not the correct ones, while the second is due to the to the additional poles introduced when 
taking the "finite-temperature" Green function instead of the vacuum one. Both of them is ultimately related to 
the zero of the function g = /32i(/32i + «7r) at P12 = itt. We would like to remark that the tanh^ factors which 
were put in by hand can be considered as arising from vertex operators of usual conformal free scalar fields with 
the two-point function — log/3i2. Since such vertex operators will have trivial braiding, one can write down a tensor 
product representation by tensoring our vertex operators with the conformal free field vertex operators. Then the 
expressions for the form-factors will automatically include the correction factors for the kinematical poles. However, 
all our attempts to find an "irreducible" (not tensor product) representation have failed so far. For this one has to 
change the g function, while retaining its property 



in such a way that we get the right operator product structure. One has to notice that the above condition determines 
g only up to multiplication by an even function. However it does not seem possible to find such a function, although 
we have no proof that such a function cannot exist. We simply have been unable to construct it so far. The solution, 
which takes tensor product representation does not seem to be superior to the procedure of simply introducing the 
additional factors. 

We have obtained a procedure with which one can (in principle) calculate form-factors of the sausage model. The 
main question, however, is to find a general procedure with which one can associate a similar free field representation 
to any factorizable S-matrix and in this way find the solutions for the form-factor axioms. The other question is 
whether the solutions obtained from the free field representation form a complete set. If this is true, then the method 
may provide a major step towards understanding the operator content of the models. We also want to mention in 
this context the relations to the work of Smirnoy HI. It would be interesting to analyse the form-factors from the 
point of KZ equations and following the lines of |10p[l2| to obtain a counting of the local fields in the model. It also 
seems worthwile to examine the rational values of ^ and to see what sort of structure emerges in that case. 



The authors would like to thank M. R. Niedermaier for useful discussion. This work was partially supported by the 
Hungarian National Science and Research Foundation. 

APPENDIX A: PROOF OF THE BOOTSTRAP FUSION AND THE ZAMOLODCHIKOV-FADDEEV 

ALGEBRA RELATIONS 

In this appendix we outline the proof of the bootstrap fusion and the ZF relation. Starting with the bootstrap 
fusion, the only nontrivial relation to prove is: 



VI. CONCLUSIONS AND FURTHER QUESTIONS 




(6.1) 



ACKNOWLEDGMENTS 
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= CZ_Z + y )Z - y ). (Al) 



Writing out in details: 



exp [gx'C/(/3) - (g + g-')xC/(/3)x + (/3)x' 



exp ( ^ (g^xV^(/3+) - g-^y(/3+)x) exp ) ((?tx^(/3-) - g-*l^(/3-)x), (A2) 



i 

Q_ 

where P± = P ± j7r/2. 

First, we prove the equahty for the integrands using the normal ordering procedure. We symmetrize in the arguments 
of the screening charges, since they are to be integrated out. We get the relation: 

iqr\f3 - li)r\P - 72) -{q + q~')9'\li " P)r\P - 72) + 9-'ff-'(7i " P)9''il2 - f3)mi2 - 71) + 
(similar term with 71 72) = 

5(72 - 7i)ti'(72 - P+)w{P^ - 7i)((7^w(/3+ - 71) - g"*w(7i - /?+))(giu;(/3_ - 72) - q~^w{-/2 - /?-)) + 
(similar term with 71 <-!■ 72) , (A3) 

which can be proven after some elementary, but lengthy calculations. Then one proves the equality the integrals, 
which amounts to checking that the contours can be deformed into each other. 

For the proof of the ZF relations one actually uses only the ratio g{l3) / g{— P) ■ Hence one can substitute the functions 
by other functions having the same ratio, eg. G in the proof. The proof is essentially the same as presented in 
In that work, we could only work out the proof up to the identities containing two screening charges. Since 
then, improving our MAPLE algorithms, we have been able to prove the ZF relations for all combinations of the Z 
operators, including those with three and four screening charges, in the general sausage model case. The procedure 
is very similar to the proof of the fusion relation, but the formulas are even more lengthy and we do not attempt to 
describe it here. For the 0(3) case it was discussed in some details in ||l[. 



APPENDIX B: EVALUATION OF TRACES USING COHERENT OSCILLATOR METHODS 

Here we show how to calculate the traces involved in the form-factors using coherent oscillator methods. We will 
take a very simple example and calculate in the 0(3) limit, when ^ = 00. We will also calculate only the non-zero 
mode part of the trace, since it is the most interesting contribution when we discuss the analytical structure in the 
main text. We will take the case when the regularizing parameter Q is finite and reproduce the function G'C) used in 
the main text. The expression we want to calculate is: 

Tr(e'^*-^e**('^i)e^'^('^=)) (Bl) 

with uj — 2tt ~ Q introduced for convenience. We omitted the indices e from K and 4> since it is clear that we use their 
regularized version and this will cause no confusion. 

We use the following redefinitions. First, we change our oscillator to have the usual normalization: 

[a„,a™] = TO(5„+„.o. (B2) 

These are related to a,n in the following way: 



a-m = V^Vm(e)am, a-m = \/ Nm{€)a-ra, (B3) 

where 

, , sinh sinh Trme t,^, 

N„,{e) = 2 e— . B4 

m 

As usual, the operators am and a_m are Hermitian conjugates. The non-zero mode part of the operator K looks 
very familiar: 

oc 

K = ie ma-mCtm + zero mode terms (B5) 

m—l 
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We normal order under the trace to get: 

5(/32i)Tr(e"*^ : e^^^^^^^e^^^*) :), (B6) 



and put in the expUcit expressions ( |3.9| ), (3.28). Observe that the trace factorizes into independent pieces over each 
harmonic oscihator. We evaluate it in the basis of coherent states. The relevant formulas are: 

a|A) = A|A), {X\a^ = (A|A, x"^"|A) = \xX), 

TrO ^ J ^e-^\x\0\X). (B7) 

Here we omitted the mode index m since the above formulas are the same for each oscillator mode. Performing the 
calculation for each mode separately, we get a Gaussian integral, which can be evaluated easily. Putting together all 
the contributions we find 

' \ sinh^ ^{1 - e--^"^) J sinh2iHil(i_e— ) j ^ ' 

Now we substitute A^m(e) and then expand 

(l_e-'^^")-i = ^(e-""^")™ (B9) 

in the exponent. We can then sum up in m using the formula 



— ^ log(l - x). (BIO) 



m— 1 



Ignoring the rapidity independent parameters (which can be assembled into the coefficient C2), we arrive at an infinite 
product formula: 

00 

J~J ^ — aen+(Tr — a)e+iemf3i2^^-^ g — Qeri+(7r — a)e— iem/3i2 ^-|^ ^ — aen—ae+ieml3i2'^^j^ g— aen— ae — iem/3i2 ^ (Bll) 

Using the definition of the deformed gamma function 

1 — exp(— 27re(a:; + n)) 

and the fact, that for e — it tends to the usual gamma function, we get for G'-<\/3i2) the expression (apart from 
constant factors): 



(/32i+«7r)sinh(^ 



r(i-!^ + !^Vfi-- + ^ 



which in the limit C ^ gives back the ^ = 00 case of (|5 



(B13) 



APPENDIX C: A SAMPLE CALCULATION OF THE OPERATOR PRODUCT SINGULARITIES 



Here we sketch the calculation of operator product singularities on the example of the product U {Pi)xxU {P2) in 
the 0(3) symmetric case. This product is given by (up to constant factors) 



(/321 + in)p2i 
1 



f dji 


f ^72 


J 2vrj 


' 2n 



(721 + «7r)72 



1 



1 



1 



(71 - /3i + «7i')(7i - Pi) (72 - l3i+ i7r)(72 - Pi) 



{P2 - 71 + »7r)(/?2 - 71) {P2 - 72 + i7r)(/32 - 72) 



(CI) 
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with the contours specified by the prescription given in the main text. We have possible pinch singularities at l32 = Pi, 
P2 = Pi — iTT and P2 = Pi — 2i7r. We can analyse their contribution by substituting Pi = P2 + ie, Pi = P2 + "tTr + ie 
and Pi = P2 + 2i7r + ie respectively into the integral and studying the limit e — > 0. One finds a regular behaviour in 
the first two cases, which is caused by the zero of the preintegral factor, cancelling the pinch contribution. Only the 
third case will give a singular contribution, which will correspond to a first order pole in the operator product. Other 
operator products can be analysed similarly. 
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